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Q ■ Abstract 

a^ : 

Q^ , This paper presents the derivation of Schwinger's gauge invariant result of 

Q I Im Ceff upto one loop approximation, for particle production in an uniform 

O^' electric field through the method of complex trajectory WKB approximation 

^ . (CWKB). The CWKB proposed by one of the author's [0 looks upon particle 

rrl production as due to the motion of a particle in complex space-time plane, 

thereby requiring tunneling paths both in space and time. Recently ||2|, ^ 
S^ . there have been some efforts to calculate the reflection and transmission co- 

H I efficients for particle production in uniform electric field that differ from our 

expressions for the same. In this paper we clarify the confusion in this regard 
and establish the correctness of CWKB. 

1 Introduction 

The complex trajectory WKB approximation has been an effective tool to under- 
stand particle production. It works in Minkowski and as well as in curved spacetime; 
both in space and time dependent gauge. Though originally proposed and devel- 
oped by one of the authors ^ , the method is now being persued by some workers 
[0, 1^, Q, technically being the same in approach with CWKB, having a slight differ- 
ences with respect to interpretation and results. To show the exactness of CWKB 
with the standard results, we consider the particle production in constant electric 
field, where gauge invariant result due to Schwinger [^ is available for comparison. 
Schwinger [^ demonstrated that the probability for a vacuum to remain a vac- 
uum in a volume L^ = V over time T is determined by the quantity exp {~2Im Ceff)VT, 
i.e., 2ImCeff is the rate of decay of the vacuum per unit 4-volume. In one loop 



approximation, the Schwinger's result for a constant and uniform electrical field is 
2Jm4V, = (2. + 1)^ E ^^e^^. (1) 



It is a gauge invariant result. Here m is the mass of the particle produced with 
charge q (for fermions {s = 1/2) the upper sign, and for bosons {s = 0) the lower 
sign is to be taken in Eq. (1)), and Eq is the electric field. The expression (1) is 
related to pair production through Nikishov's representation : 



2Jm£«L3T = TE / 7^^'ln(l±^p), (2) 



(2^ 

Up = exp (-TT — ^), (3) 

QEq 

where Up is the average number of pairs formed by the field in the state with a given 
momenta and projection of spin k = k,s. 

The above technique of obtaining the Im Ceff can be generalized to arbitrary 
electromagnetic and gravitational backgrounds but is an uphill task since it is not at 
all convenient to regularize the effective Lagrangian in all such cases. Due to these 
reasons many workers in this field turn towards the Bogolubov transformation tech- 
nique based on the method of normal mode analysis to understand the phenomenon 
of particle production. Here the problem arises with respect to interpretation when 
one tries to obtain the gauge invariant results, mentioned in (1) and (2). In time 
dependent gauge the instability of the vacuum is quite understandable. The vacuum 
|0_ > at t — > — cxD is not the same as the vacuum |0+ > at t — > +cxd. The vacuum 
decay |0„ >^ |0+ > allows one to see particle production in which the Bogolubov 
coefficient Px (A is a parameter that specifies the state of produced particles) is not 
zero and is interpreted as 

<0_|iV(A)|0+>=El/5Ar, (4) 

A 

where N{\) is the number of 'in' Wj-mode particles found in the out vacuum |0+ >, 
where Ui{\) modes construct the in- vacuum |0_ > and 'Uj(A) modes construct the 
out- vacuum |0+ >. When one considers the same problem in space dependent gauge; 
the vacuua at a; — ^ =Fcxo remain the same and there is no vacuum instability. Con- 
sequently one expects no particle production in a space dependent gauge. However 
through mode analysis one still recovers a non-zero |/3aP even in space dependent 
gauge, implying particle production. A tunneling interpretation is then invoked to 
interpret particle production in space dependent gauge because of non-zero |/?aP- 
In this interpretation one switches to the interpretation of 'tunneling through the 
barrier' [Q to calculate the reflection and transmission amplitude R and T and to 
obtain |/?;^p through a prescription, identifying |aAp =t \Px\^ = 1 as being equivalent 



to |Tp =F |/2p = 1 (± refer to spinor/scalar fields). The prescription that relates 
R and T with a\ and (3\ is found to work differently in space and time dependent 
gauges. Recently a complex path analysis based on a technique of Landau [0] is 
being used 0, || to obtain (3x from the expressions of R and T. In CWKB, we 
used the technique of ref. [§ using the idea of Cornwall and Tiktopoulos 0. As 
pointed by us earlier, this has now led many to realize the usefulness of complex 
path analysis to understand particle production. In our view, the analytic contin- 
uation inherent in Bogolubov's technique is exploited in ref. [Q with an advantage 
towards the point of interpretation. However these works do not clarify the source 
of instability, how and where the particle production occurs. In these works 0, ^ ^ 
one also requires an advance knowledge of the gauge invariant result to settle the 
relation between the pairs (i?, T) and {a\,l3\). The knowledge of mode solutions 
for the analytic continuation x -^ — oo to x — * +oo (or, t -^ — cxo to t — > +oo) is 
a basic requirement in these works. In principle one should obtain A^(A), without 
connecting it to I3\, from the results of R and T themselves. There is another point 
of significant argument that has motivated the present work. In CWKB, the results 
of R and T do not coincide with ref. 0, ^, Q and hence needs discussion. In ref. 
[0, |], attempt has been made to obtain the results (1) and (2), but being partially 
successful. In the present paper we use CWKB to obtain the reflection amplitude 
-RcwKB and transmission amplitude Tcwkb and use a first principle derivation of 
N{X) without reverting to the calculation of Bogolubov co-efiicient and obtain the 
exact results, term by term, of (1) and (2) along with the prefactor in the expres- 
sion (1). In CWKB, we do not require to know the mode solutions at t — ^ ±oo 
(or X — » ±oo). It is observed that -Rcwkb is related to the full S"- matrix elements 
whereas Re, in the works of ref. [Q is related to the connected 5*- matrix element. 

In CWKB, the particle production is considered as a 'reflection' process either 
in time (in time dependent gauge) or in space (in space dependent gauge). This is 
the well known Klein paradox like situation, found in all standard textbooks where 
the pair creation occurs in a barrier region due to reflection causing Zitterbewegung 
within wavepacket in the barrier region. Only when the potential is capable of 
transferring energy at least 2mc^ to the wave packet, i.e., the minimum energy 
required to conserve energy in making a pair, will then real "permanent" particles 
be created with the packet moving past the potential. The Zitterbewegung region is 
specified by two points within which the mixing of positive and negative frequency 
parts takes place and the interference term varies in time with a frequency at least 
2mc^. The 'two points' are the turning points from where the particle turns back. 
The turning points for a given problem is fixed by the one dimensional Schroedinger- 
like equation either in space or in time depending on the gauge chosen. 

The present work is motivated to deal with particle production through the 
technique of complex trajectory WKB (CWKB) approximation, proposed by us in 
various earlier works [|TU], |TT], |T2|, [T^ (mainly restricted to particle production in 
expanding spacetime). The organization of the paper is as follows. In section II 
we discuss the particle production as a process of reflection in time. It is a Klein- 



Paradox-like situation not in space but in time. We choose an idealised example to 
elucidate our approach. In section III we show a method to calculate the number of 
produced particles in terms of reflection coefficient and obtain the Schwinger's result. 
In section IV we obtain the relation connecting the Bogolubov coefficients and the 
reflection coefficients. The concluding section synopsizes the results obtained in the 
present work and clarify some misconception on the CWKB results. 

2 Basic Principles of CWKB 

The method of CWKB in the context of cosmology was proposed by one of the 
authors in [|I[ and subsequently developed in [^, |n|, |r^, [TB|, ^^ [T^ to treat particle 



production in Robertson- Walker spacetime and deSitter spacetime. In Feynman 
spacetime diagram a pair production is looked upon as Vac —>■ e~^ + e~, in which the 
vacuum supplies energy to create the pair 2mc^ such that the positron e+ of energy 
Ee+ and the electron of energy Eg- both move forward in time. Using Feynman- 
Stuckleberg prescription, the electron moving forward in time can be considered 
as positron of energy —E^- moving backward in time. The net result is that the 
pair production is looked upon as a process of reflection, not in space but in time 
from vacuum 0. The change —E^- to +Ee- requires that the system showing pair 
production through process of reflection must have turning points in time and as 
the system evolves from t -^ — oo to t ^ +00, there is rotation of currents showing 
antiparticle to particle rotation. Obviously the turning points are complex in time, 
reflecting the situation of 'over the barrier reflection'. 

To get an idea of the above mechanism consider an one dimensional Schroedinger 
like equation in time 



Eq.(5) has turning points given by 



+ u;\t)ij = 0. (5) 



a;'(ti,2) = (6) 

where the turning points ti^2 are complex. Balian and Bloch [^ formulated an 



approach to quantum mechanics starting from classical trajectories with complex 
coordinates. In the approach, within the framework of WKB approximation, wave 
optics is generalized to complex trajectories to build up the wave. Such a method 
accounts for contribution of exp (— f ) to the usual WKB wave and is found to repro- 
duce quantitively all quantum mechanical effects even in cases where the scattering 
potential varies rapidly over a distance of wavelength or less. We applied this tech- 
nique following [§], to obtain an expression for ip(t) at t — >■ co as 

i[j(t) ~^ — L^ Usit,to) _ ,jie-iS(t,to)] (7) 



where 

S{t,to) = [ uj{t)dt (8) 

Jto 

and 

^ exp[2iS{ti,to)] 

l + exp[2tS{h,t2)y ^ ^ 

The detailed derivation of Eq. (7) will be found in |^ and also in our previous 



works [jT2|, [T3[. The interpretation of Eq. (7) is as follows. A wave starts at to > ^ 
and moving leftward arrives at t. This is represented by the first term in Eq. (7). 
Another wave starts from Iq, moving leftward arrives at the turning point ti to get 
reflected to arrive back at t. This is a reflected wave contribution given by 

exp [+iS{ti, to) - iS{t, ti)] = exp [2iS{t, to) - iS{t, to)]. (10) 

The contribution (10) is then multiplied by repeated reflections between the complex 
turning points ti and ^2 and its contribution is 

OO 1 

The contribution (10) times (11) give the second term in Eq. (7) as the reflected 
wave where R is the reflection amplitude. The essence of Eq. (7) is that there is no 
particle at t ^ — oo i.e., 

^|J^n^^^exp[^Sit,to)], t<0 (12) 

but at t ^ +00, Eq. (12) evolves into 

^i^ t ^oo (c^(t))i/2 [^^P t^*^^^' ^o) ~ iR[-iSit, to)]]] . (13) 

We identify R as pair production amplitude. The above method works also in space 
dependent gauge in which we are to replace Eq. (12) and Eq. (13) as 

X -^ oo ^ ' 

X -^ — oo 

We applied this technique for particle production by external electromagnetic back- 
ground for A^ = {EqX, 0, 0, 0) as follows. The Klein-Gordon equation is 

{{df, + iqA^){d^ + iqA^')+m^)^ = 0. (16) 

With M = kl + k"^ + m^ = kj_ + m"^ and 



we find from Eq. (16) 



d' 



2, 



dx^ 



+ 



u + qEof-M^ = 0. (18) 



Let us substitute 



LJ 



in Eq. (18) to get 



dh 



p = ^qEoX + -j=, (19) 

kl + kl + m^ 
A = ^--^ , (20) 

Q^2 + (P' - ^)^ = 0' (21) 

so tliat tlie turning points are at p = ±-\/A. We evaluate R using Eq. (9) and Eq. 
(20) as 

We evaluated tlie same problem in time dependent gauge, found complex turning 
points in complex time plane to obtain the same expression of R. Thus the reflection 
amplitude R in CWKB is a gauge independent result. 

Recently |0, ^ a complex path analysis similar to us [they used the technique of 
Landau and needs to know the mode solutions] obtained the result 

„-7rA 

l^cl' = — r- (23) 

This result is also obtained in 0] but differs from our expression (22). Henceforth 
we call our reflection amplitude as -Rcwkb to distinguish it from Eq. (23). In 0], 
the unitarity relation (that also follows from charge conservation) 

|i?c|' + |Tcr = l, (24) 

is equated with 

\ax\^ - \Px\^ = 1 (25) 

(here ax and Px are Bogolubov coefficients) to obtain 

i/^^r = if (26) 



\Tc 
1 



l«Ap = ^. (27) 

No arguments are placed to obtain Eq. (26) and Eq. (27) though |/3aP = |^ 
1 



IRcP 
and I^Ap = tr^ is also a possible choice. It is not at all clear what transpires to 



connect |-RcP and |Tcp with \a\\^ and |/?aP, wherein one can in principle determine 
the number of pairs formed (= \I3\\^) simply from the results of |i?p only. The 
reason for such a choice may be that one knows a priori the result of Schwinger to 



settle Eqs. (26) and (27). In a recent work |T^ we call into question the validity 
of Eqs. (26) and (27). In the present work we show that both \R\ and |i?c| given 
in Eqs. (22) and (23) are correct, in which i?(= full S-matrix element) is related 
to disconnected propagator and Rc{= connected S-matrix element) is related to 
connected propagator. We also derive the relations (26) and (27) from our result of 
R. 



3 Reflection Coefficient And Average Number Of 
Pairs 

The amplitude for observing a scalar particle at the point Xa and an antiparticle at 



Xb is given by [18, 19 



A = -Aq / dkadkbda^da^*{ka, Xa) d^a G{xb, Xa) d^b f*{h, Xb) (28) 

where da^ is the element of a space-like hypersurface, 9^ is the derivative in the 
time like direction and Aq is the amplitude for no particle production. Consider 
scalar particle in time dependent gauge. In such cases, the x dependence of / and G 
are all phase factors and hence the integrations over d^Xa^b (required in evaluating 
the Green function G{xb,Xa)) and d^ka^ fix the above numbers of the particle and 
antiparticle to be —ka = kb = k. If we look at ta = tb, we see a particle and an 
antiparticle on the same space-like surface with opposite momenta. Under these 
circumstances Eq. (28) can be written as 

A=fd^kA{k) (29) 

where A{k) being the amplitude for creation of a pair with momenta k and — k. For 
an electromagnetic potential A^ = (0, ii^o^^O, 0), the Klein-Gordon equation is 



m'f = i-d^ + (d^ - qEoty + d^ + dl)f. (30) 

1 
/(x,t) = -- 

we find 



y 

With 

1 

Ak.x 



fi^^t) = -^-^^e^'^-g{t), (31) 



9 = ND_^_,.{^2z^o{t - ^), (32) 

where 

\N\ = (2gEo)-^e-'^, (33) 



and 



\^^{m' + kl + kl). (34) 



Evaluating the propagator it can be shown that A{k) turns out to the form 

\A{k)\' = \Ao\'uk (35) 

with 

'^^ = TT^- (36) 

Thus |^(/c)| is the probabihty for the creation of one pair with wave number k, |Aop 
is the absolute probability of creating no particles in a given mode (i.e., in mode 
k) and Uk is the relative probability of producing a pair if none were present in the 
initial state. This ujf^ is the square of the connected S-matrix element ||^ and is 
denoted by -Ke- 
lt should be pointed out that in the complex path analysis, we could start from 
the relation (35) itself to determine |v4op and u^ using the result of -Rcwkb with 
|y4(A;)p = |i?cwKBp- Henceforth we omit the subscript 'CWKB' i.e., write -Rcwkb = 
R. It may be pointed out that R is the full ^-matrix elements whereas the Uk by 
definition is related only to the connected ^-matrix elements so that uj^ = \Rc{k)\'^ 
i.e.. Re is related to the connected propagator Sac whereas R is related to the 
disconnected propagator Sa- These propagators are related by 

Sa = e^^^Ac (37) 

where W = Tr IuS'ac is the sum of connected vacuum graphs [see 0]. 

The probability of producing n pairs of particles with wave number k is then 

P„(fc) = \A,\'u;l 

Hence the total probability of creating n = 0, 1, 2, . . . pair will sum to unity (for 
Bose statistics) i.e., 

\AoWl + u,^ + ul + ...) = 1 

or, \Ao\^ = 1 - ^fc- (38) 

Hence the absolute probability for 1-pair production is 

\A{k)\' ^ P,{k) = {1 - u^P^,. (39) 

By definition the reflection coefficients |-RcwkbP = l-Rp equals Eq. (39) so that 

\R\^ = Uk{l-uJk). (40) 



The average number of pairs with wavenumber k is 

N{k) = J2nPM = -r^^- (41) 



Let us now relate Eqs. (38), (40) and (41) to the vacuum persistence probabihty. By 
definition \Aq\'^ is the probabihty of no pair production for wave number k. Hence 
the total probability of no pair production for all wave number k is 



\-i 



= e-Efei°(i+^W). (42) 

Eq. (42) gives the probability that no particles are created in the vacuum in any 
given mode i.e., the probability that the vacuum remains the vacuum is given by 

I < out\in > P = J||Ao(A?)P 

k 

or, e-2^™^=^^-^^ = e-Sfc'''(^+^(^)) 

or,2ImCeffV4 = J2^n{l + N{k)), (43) 

k 

where 14 = L^T. In CWKB 

-nX 

[l + e 
Using Eq. (40), we find two roots 



,2 e 



Rr = Tr--^r (44) 



1 e-^^ 

1 + e-'^^ ' 1 + e- 



^i = r, r- (45) 



We take the second root so that 

„-7rA 

l + e 



^t: = r- (46) 



l-OJr= i -. (47) 

The reasons for such a choice will shortly be clear. Using the values of Eqs. (46) 
and (47) in Eq. (41), we get 

N{k) = e-^^ (48) 

so that 

2ImCeff.Vi = ^ In (1 + e""^). (49) 

k 

9 



If the electric field Eq -^ 0, \out >^ \in > and there is no instability so that 
< out\in >= 1. This implies that ImCe/f — * 0. From (34), we see A = ^ ^ -^ cxd 
when Eq -^ 0. Hence from Eq. (40) we have ImCeff = implying < out\in >= 1. 
Had we chosen the other root in Eq. (45), we would get N{k) = exp (ttA) and it 
goes to infinity as Eq — > 0, giving < out\in >^ infinite. This is unphysical. Thus 
the correct choice is given by Eq. (46). 

Expanding the log term in Eq. (43) we get 

V,.2ImCeff = 5](-l)«+iie-'^"^ 

k,n 

= y(_i)«+iie — ^0 — . 
"r^ n 

k,n 

We convert the sum by integral as follows 

- ^-%^U-^"*'^^^- (50) 

Since V4 = L^T, we get 

2/™A„=(^|(-ir«-ie^^ (51) 

This is exactly the Schwinger result for ImCeff in one loop approximation. Eq. 
(43) is the Nikishov's representation (Eq. (2)). Our derivation no where requires 
to evaluate the Bogolubov co-efficient or to know the mode solutions. This amply 
refiects the gauge invariant content of CWKB. 

4 Reflection Coefficient And Bogolubov Coeffi- 
cients 

Let the positive and negative frequency solutions be denoted by 0fc± (in the limit 

t -^ —00) and 

k 

af\Vac>=af\in>=0. (53) 

Eq. (53) specifies the in vacuum. For out vacuum let us define 

'^=E(«r4^^+^r4"^ (54) 



10 



af^out >= (55) 

where 0^ are the positive and negative frequency solutions at t ^ +oo. The 
Bogolubov coefficients a^ /9g are now given by 

•^M^) = "^""^ + /^4"^ (56) 

*^H = /5|</>^^ + «l<^lr^ (57) 

with 

i^d' - m' = 1, 
4"^ = [4"^]*' <^k.) = i<^H-)r- 

The operators a^ and 6^ are connected to the out operators as (we now omit the 
vector sign over k) 

aT = alar' - PlhT\ (58) 

6+- = -/3,ar* + afc^r"*- (59) 

The pair formation amphtude is 

apair =< out\ak'outb^out\'^^ > ■ (60) 

Using (59), we write it as 

3* 
apair = < out\ak'outialr\-^ + ^ato^t)\in > 

3* 
= -^ < out\ak'outatout\^i^ > 

O-k 

8* 8* 

= -^ < out\6kk' + a'^outO'k'outlin >= -7 < out\in >, (61) 

since [akoutydt' out] — ^kk' and < out\al = 0. The relative amplitude of producing a 
pair if none were present in the initial state is 

< out\ak'outbkout\in > , . 

< out\in > 

The square of this amplitude gives the probability of creating a pair in the given 
mode. It is Uk as denoted earlier [see Eq. (36) and the discussion thereafter]. Hence, 
using (61) and (62) 

\Pk :2 I D |2 



cuk = \^\'=\Rc\. (63) 



We thus get 



UJk \Pi 



'k\ 



|2 



l-UJk IftfcP - |/9fc| 

11 



\Pk\'. (64) 



Hence 



e 



-ttA 



Ificl^ = YT^. (65) 

This result of Rq was derived in 0, |] using the technique of Landau and also in 
[§] by saddle point method. Let us try to understand the relation between R and 
Re. 



|R|2 — IP l2 

\K\ — T- -rtf 



1 

"tcl 



^ ln(l+iV{fc))|D |2 

= e-''^"'^^fi<^^-^'\Rc\^. (66) 

This connection between R and i?c was also cited in ref. ||^ where 

2Im CeffVi = 2y^^ Im Cef f{k) .Vj 

k 

= ^ln(l + Ar(A;)). (67) 



From (63), we find 



l^cl^ = 1^, (68) 



so that 

|2 1 I D |2 l^fef ~ \Pk 



cl =l-|^c| = f7-[^ =1 — 1^- (69) 



\Oik\ l«fc 

From the unitarity relation |-RcP + l^cP = 1; we get 



1 



Hence 



\Tcr = r^- (70) 

I D 12 1 

1/312 r^-C I |2 -■- /-TI \ 

The results (70) and (71) coincide with (26) and (27) as were taken in ^. The 
CWKB result of R thus correctly reproduces the result of Rq obtained by other 
methods. 

5 Conclusion 

Let us try to synopsize the results. In the complex path analysis [0, that uses 
Landau's technique or in the saddle point method Q, one obtains Rq and Tc- 
These Rq and Tc are related by charge conservation as follow 

|i?cP + l^cP = 1; {scalar case) (72) 

12 



1 + l-RcP = l^cP, {spinor case). (73) 

Using the unitarity constraint of Bogolubov coefficients, one determines 1/3^^ with 
the identification N{k) = |/3fcP- Using Nikishov representation one relates N{k) with 
ImCeff to obtain the Schwinger resuh. Here one needs to know the mode solutions 
as well as, the continuation from t -^ — oo (x -^ — oo) to t ^ +oo (a; -^ +oo). 
Moreover in 0, 0], there is an ambiguity in fixing (3k and ak from Eqs. (72) and (73). 
In CWKB, we calculate R by simple integrals and particle and antiparticle states are 
fixed by WKB definition (Parker's definition). This fixed the vacuum. Nikishov's 
and Schwinger's results follow through our derivation. The instability of the vacuum 
that causes particle production occurs due to motion in complex-coordinate plane. 
There is a rotation of currents as the system evolves from t -^ — oo [x —>■ — oo) 
to t ^ +00 {x — * +oo). The mixing of positive and negative frequency states 
occurs between the turning points. The value of the wavefunction at real t is also 
contributed, not only by real trajectories but also by complex trajectories, in the 
sense of WKB approximation. Our result shows that CWKB justifiably takes at 
least one loop quantum correction within the semiclassical approximation. We have 
shown elsewhere in detail that the particle production in CWKB is equivalent to 
tunneling from left Rindler wedge to right Rindler wedge via real x and imaginary 
time plane. Not realizing the difference between R and Re and the beauty of the 
complex path analysis, some circles argue the expression of R as wrong comparing 
it with Re- We clarify our view in this respect and establish the correctness of 
CWKB. 
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